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Iron abundances from optical Fe III absorption lines in 
B-type stellar spectra 
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ABSTRACT 

The role of optical Fe ill absorption lines in B-type stars as iron abundance diag- 
nostics is considered. To date, ultraviolet Fe lines have been widely used in B-type 
stars, although line blending can severely hinder their diagnostic power. Using optical 
spectra, covering a wavelength range ~ 3560 - 9200 A, a sample of Galactic B-type 
main-sequence and supergiant stars of spectral types BO. 5 to B7 are investigated. A 
comparison of the observed Feiii spectra of supergiants, and those predicted from 
the model atmosphere codes tlusty (plane-parallel, non-LTE), with spectra gener- 
ated using SYNSPEC (LTE), and CMFGEN (spherical, non-LTE), reveal that non-LTE 
effects appear small. In addition, a sample of main-sequence and supergiant objects, 
observed with FEROS, reveal LTE abundance estimates consistent with the Galactic 
environment and previous optical studies. Based on the present study, we list a num- 
ber of Fe III transitions which we recommend for estimating the iron abundance from 
early B-type stellar spectra. 

Key words: atomic data ~ line: identification - Galaxy: abundances - stars: abun- 
dances - (stars:) supergiants - stars: variables: other 



1 INTRODUCTION 



Iron lines dominate the 
phys ical objects, such 
I1997I : iHatzidimitriou et all 
II reRions IIRubin et al 



spectra 



of man y astro- 
(jMcKenna et al.l 

photoioni zed H 

iRodri'guej |2002| : 

lEsteban et al. 2002 ) and active galactic nuclei 



JSigut fc PradhanI l2003l: ISigut. Pradhan fc Nahaj |2004 



IZhang. Dultzin-Hacvan fc Wand 120071 ). The atomic pro- 
cesses for iron and other iron-group ions have been the 
subject of numerous investiga tions, for example the IRON 
Project (|Hummer et al.|[l993l ) which consid ers applications 
in as trophysical and laboratory plasmas (jPradhan et al.l 
Il996h . Absorption lines of iron provide important metal- 
licity diagnostics for both stars and galaxies, and also 
play a key role in investigating star fo rmation histories 
throu gh element ratios, such as [a/Fe] j Cilmore fc W ys3 
Il998( ). However, there is a lack of robust Fe abundance 
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torio del Roque de los Muchachos of the Instituto de Astrofisica 
de Canarias 



determinations in exte rnal galaxies, e.g. the Magellanic 
Clouds (see for example. jRolleston. Trundle fc Dufton,.200^ : 
iTrundle et al]|2002l . |2007| and iMokiem et al.ll2007^ ■ due to 
the complex Grotian diagrams for Fe, the low metallicity 
environment of the Magellanic Clouds and the reliability of 
the currently available atomic data. 

Early B-type stars are important for study- 
ing the chemica l com p osition of ou r own and other 
galaxies (iKilianI Il992l : iDuftonl Il998l ). In the optical 
spectra of B-type stars, iron lines due to a number 
of ionization sta ges are observed (see for example, 



Gies fc Lambert] 

I993I : ISmartt. Dufton fc Lennon i997i; iMorel et al.l 



19921: II ennon, Dufton fc Fitzsimmond 



2006 



Absorp tion features arising fro m Felll are primarily de- 
tected (|Hardorp fc S cholj Wol), with Fell also found in 
later B-type stars (Pintado fc AdelmanI Il993l ). Feiv lines 
are not expected in the optical spectra of B-type stars due 
to their intrinsic weakness in this temperature range. 

The optical Felll line spectrum has not been widely 
employed to determine abundances, as it has often been 
believe d to be too weak to provide reliable measure- 
ments l|Kendall et al.l Il994al ). However, it has been used 
for chemical composition studies of several bright, narrow- 
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lined main-sequence B -type stars, such as C Gas, 7 Peg , 
tHer, rSco and ALep dSniiderdllQegl: iPeters fc Allej|l97Cll: 
Hardorp fc Schold Il970l : IPeterj Il976l : IPeters fc Polidanl 



19851 '). These Galactic objects have sufficiently high Fe 



content, coupled with narrow metal absorption lines (due 
to their low projected rotational velocities), so that, 
even with relatively poor quality optical spectra, Felll 
features can be detected. However, the quality of the 
available spectra did not allow the study of Felll lines 
in objects of lower metallicity. Therefore, more recently 
the optical wavelength range has been largel y over- 
looke d in favour of the u lt raviol et domain ( Swings fc Vreuxl 
1 19761 : iPeters fc Polidanl Il985l : iDixon fc Hurwit j Il998h 
and in pa rticular the very rich spectral r egion around 
1900 A (iThompson. Humphries fc Nandvl 1 1974 iHebeij 



iiJUU A (iinompson. numpnries iNandvi iMeber 
19831 : 'Kenda ll et al l Il994al : brigsbv. MuUiss fc Baeij 1 19961 : 



Moeh lcr ct al. 19981 ). On the other hand, due to the high 



density of absorption features and resultant blending, con- 
tinuum placement in the ultraviolet is difficult, and hence 
significant errors ma y be present in the derived abundances 
l|Moehler et al.lll998l ). 



Due to instrumental advances in more recent years, it 
has become routine to obtain high resolution and signal-to- 
noise (S/N) spectra, and the use of the optical Felll lines 
as a diagnostic ha s been revisited (e.g. iGunha fc Lambert) 
1 1994 iKihanI Il994l ). A number of main-sequence objects 
have subs equently been re-examined, for example 7 Peg 
andt Her (jPintado fc Adelmanlll993l : E iu. Sterken and Zhao! 
1 19981 ). producing results consistent with the earlier optical 
analyses. In addition, the higher S/N ratios of the available 
spectra have made it possible to detect the weak optical 
Felll lines in a number of lower metallicit y objects, such 
as AV 304 in the Small Magellanic Gloud ( RoUeston et al.l 
I2OO3I '). an d the globula r cluster post- AGE stars ZNG-1 
in MIO (|Moonev etHI [200 4). J3arnard 29 in M 13 and 
ROA5701 in tjGen ijThompson et al.|[20d7l '). 



A number of B-type stars have been studied using 
both ultraviolet and optical spectra, but the correspond- 
ing abundance estimates are generally in poor agreement. 
Estimates from the ultraviolet spectra are consistently 
lower than those from the optical, for example in the 
post- AGB stars Barnard 29, ROA5701 llThompson et all 
I2OO7I ). BP +33°2642 (Na piwotzki. Heber fc KoppenI 11991 ) 
and HD 177566 (iKendall et al.l ll994a"bV You ng Galactic 
objects includ i ng 7 Peg i Moehler et al., Il998l ) and tHer 
ijGrigsbv et al.l Il996l ) , plus stars in known metallicity en- 
viron ments such as the Magellanic Glouds (jPufton et al.l 
bOOTl ). also display similar discrepancies. 



Here, high resolution spectra for a number of narrow- 
lined Galactic B-type main-sequence and supergiant stars, 
covering a range of spectral types, are analysed using LTE 
and non-LTE model atmosphere techniques. Details on the 
observations, models and Felll line selection can be found 
in Sections [2] and [3] An LTE approximation is considered, 
and the reliability of such an assumption, along with an as- 
sessment of the individual Fe iii lines, is discussed in Section 



2 OBSERVATIONS AND DATA REDUCTION 

The observational data included here are used for two inves- 
tigations. Firs t ly, to allow a comparison between the tlusty 
(lHubenvll988l : lHub eny fc Lanz 1995:' Hubenv. Heap fc Land 
I199I) and CMFGEN (Hillicr ct al. 200;^ model atmosphere 
codes, and secondly, an investigation of the reliability 
of the Felll spectrum to estimate abundances. The first 
investigation utilizes a subse t of t he stars discussed in 
ICrowther. Lennon fc WalbornI (|2006') for which there were 
suitable data from ISIS /WHT and FEROS/ESO, and these 
objects are listed in Table [1] ICrowther et all l|2006l ) anal- 
ysed a sample of 25 Galactic B-type supergiants using op- 
tical spectra from three sources; CTIO data obtained us- 
ing the Loral 1 K spectrograph, WHT spectra taken us- 
ing the ISIS spectrograph (Smartt, private communication) 
and JKT spectra from the Richardson-B realey spectrograph 
(jLennon. Dufton fc FitzsimmonsI Il992l ). A number of ob- 
jects were omitted from this study due to having poorly 
detected or unseen Felll features (i.e. objects of spectral 
type BO. 5 and earlier), or if their atmospheric parameters 
lay outside the available tlusty grid points (for example, 
HD 13854, HD 152236, HD 190603 and HD 194279). For the 
second investigation, a sample of B-type stars, consisting 
of two main-sequence stars and nine supergiants, has been 
selected, from a FEROS/ESO dataset. These are analysed 
using the tlusty atmosphere code and are listed in Table 
[2| The objects have been chosen to cover a range of spectral 
types and have relatively narrow-lined absorption spectra. 
The supergiant HD 53138 has been included in both samples 
due to it having a rich Felll line spectrum. 

The WHT data were taken using the Intermediate- 
dispersion Spectrograph and Imaging System (ISIS), on 
1999 July 08. The EEV12 CGD detector (2148x4200 pbcels) 
was used on the blue arm with the R1200 grating and a slit 
width of 1.0 arcsec. This resulted in a wavelength coverage of 
~ 3920-4720 A. Typical S/N values for spectra at ~ 4500 A 
are given m Table [U The data were reduced using standard 
procedures within iraf, which included flat-fielding, bias 
subtraction and wavelength calibration and the combined 
spect ra read into the Starlink package dipso (jHowarth et al.l 
I2OO4I ) for further analysis. 

High resolution ESO spectra were obtained using the 
Fiber fed Extende d Range Optical Spectrograph (FEROS, 
iKaufer et al.ll 19991 ). between 2005 AprU 21 and 24 (proposal 
075.D-0103(A)). The EV 2Kx4K CGD detector and 79 lines 
mm~^ echelle grating were used, resulting in a spectral res- 
olution of i? ^ 48000 and a wavelength coverage of ^ 3560- 
9200 A. Data were reduced using the standard FEROS Data 
Reduction System (DRS), implemented with ESO-midas 
software. The reduced spectra for a given target were sub- 
sequently combined using SCOMBINE within iRAlQ, and im- 
ported into DIPSO for further analysis. Their average S/N 
ratios at ~ 4500 A are given in Tables [1] and (2] 

The spectra were normalised with low-order polynomi- 
als and radial velocity shifted to the rest frame using un- 
blended metal and hydrogen lines. Absorption line equiv- 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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20071 ). 



alent widths were estimated by fitting Gaussian profiles 
to the lines. Tlie er ror in measuring these is discussed in 
iDufton et af] l|l990l l. with a well-observed unblended fea- 
ture accurate to 10 per cent (designated 'a' in Tables |4l [5] 
and [6]), weak or blended features accurate to 20 per cent 
(designated 'b') and weaker features accurate to less than 
20 per cent (designated 'c'). 



3 DATA ANALYSIS 

3.1 Model atmosphere techniques 

Two model atmosphere codes have been used in this paper, 
namely tlusty and cmfgen, which are discussed below. 

Non-LTE model atmosphere grids, calculated with 
TLUSTY and SYNSPE C (|Hubenvlll988l : iHubenv fc Lanzlll995l : 
iHubenv et al] Il998l ). have been used to derive atmo- 
spheric paramet ers. The me t hods i nvolved are d is cusse d 
in, for example, iRyans et al.| l|2003l ) . [Hunter et alj 
[Thompson et al.l (2006, 2003) and iTrundle et Zf 
wh ilst greater d etail about the grids may be found 
in iDufton et al.l l|2005l ) and at http:/ /star. pst.q ub.ac.uk, 
Briefly, four grids were calculated corresponding to the 
Galaxy, the Large (LMC) and Small (SMC) Magellanic 
Clouds and a lower metallicity regime. The Galactic grid has 
been adopted throughout this analysis and covered a range 
of effective temperatures (12000 to 35000 K), logarithmic 
gravities (log g — 4.5 dex down to the Eddington limit; g 
in cms~^) and microturbulences (^ = 0, 5, 10, 20 and 30 
kms~^). Typical spacing between grid points for the effec- 
tive temperature and surface gravities were 2500 K and 0.25 
dex, respectively. For accurate profile fitting of the hydro- 
gen lines, further models were executed at a resolution of 0.1 
dex for log g, with these models used only in determining the 
surface gravity. The metal line blanketing was assumed to be 
dominated by iron, with the exclusion of nickel in our models 
unlikely to be a source of significant err or in the atmospheri c 
parameters and chemical composition (|Hubenv et al.lll998h . 

Iron is included in the non-LTE tlusty model atmo- 
sphere calculations with sets of ionic levels being collected 
into superlevels in order to make the calculations tractable. 
Although in the spectral synthesis code it is possible to as- 
sign ionic levels to appropriate superlevels, we have instead 
opted to calculate the Fe ill spectrum in an LTE approxima- 
tion. These were then employed with the measured equiva- 
lent widths to deri ve the abundance estim ates shown in Ta- 
bles m [5] and [6] (see [Thompson et al.l 2007 ). Atomic d ata fo r 
the Felll lines were taken from iNahar fc PradhanI (|l996l ). 
as these are (to our knowledge) the most up-to-date val- 
ues available. Their oscillator strengths are shown in Ta- 
ble [3] (available in full on-line), along with t hose from the 
Kurucz d atabase (http: /7nova. astro. umd . edu , ^ubenvl 19881 : 
lHubenv~fc Lanz 1995; Hubcny ot al. 1998), plus the wave- 
length, transition, and g; and g_, values for the selected iron 
lines. 

The TLUSTY model atmosp here grids are based upo n 
the standard solar composition (|Grevesse fc Sauval|[l99i ). 
while SYNSPEC allowed a range of iron abundances to be 
calculated. Thus, using the Galactic grid, spectra were syn- 
thesi sed for a ran g e of i ron abundances from 6.7 and 8.3 
dex. iDufton et all l|2007h discuss the possible inconsisten- 
cies which can arise from the adopted abundances used in 



both the model atmosphere and spectrum synthesis calcu- 
lations, i.e. the use of SYNPLOT to synthesize spectra for 
a range of metallicities which differ from that of the un- 
derlying model atm osphere calculation. Tests discussed by 
IDufton et al.l (|2007| ) indicate that these should not be a sig- 
nificant source of error, and the method has been widely used 
(see for e xample KilianI 19941 : iMcErlean. Lennon fc DuftonI 
Il999l and lKorn et al.ll2000l '). 

CMFGEN is a non-LTE line blanketed multi-purpose at- 
mospheric code designed for the spectral analysis of stars 
with stellar winds. It solves the radiative transfer equa- 
tion for spherical geometry in the co-moving frame, un- 
der the constraints of statistical and radiative equilibrium 
iHillier fc Milleij|l998l ). As CMFGEN does not solve the mo- 
mentum equation, the mass-loss rate and velocity law must 

be s pecified. 

ICrowther et al.l ||2006|) used the current version of CMF- 



. . 20031) and a gain adopted the atomic data 

20051) . of iNahar fc PradhanI p"99il_ for the Felll lines. A model 



GE N jHillier et al 



similar to that applie d by Evans et al.l l|2004h was used by 
ICrowther et all l|2006l ). comprising 26 ions of H, He, C, N, O, 
Mg, Al, Si, S, Ca and Fe for early B subtypes, containing a 
total of 2636 individual levels which were grouped into 766 
superlevels, with a full array of 25960 bound-bound tran- 
sitions. Higher ionization stages were excluded for spectral 
types B2.5 and B3. The temperature structure was deter- 
mined by radiative equilibrium, with the velocity parameter- 
ized with a classical /3-type law for the supersonic part (with 
the exponent in the range (3 — 1-3, selected on the basis of 
the Ha profile), and the subsonic velocity set by the corre- 
sponding H-He TLUSTY photospheric model s calculated for 
each CMFGEN model atmosphere calculation (|Hillier fc Land 
120011 ). Test calculations have been performed for a H-He 
TLUSTY photospheric model and a blanketed tlusty model 
(Dr. S. Searle, private communication). The results showed 
comparable stellar parameters, and hence the use of the H- 
He TLUSTY model here appears valid for this analysis. Syn- 
thetic spectra for each star were obtained for iron abun- 
dances of Fe = 7.35, 7.65 and 7.95 dex, adopting the appro- 
priate atmospheric parameters. These were used to derive 
the abundance estimates given in Table |3] (see Section 13.31 
for further details). 

3.2 Atmospheric parameters 

The atmospheric parameters of the objects used in the com- 
parison between the model atmosphere codes CMFGEN and 
TLUSTY are listed in Table [T] Table [5] contains the atmo- 
spheric parameters of the additional objects which have been 
ana lysed solely with tlusty. 

ICrowther et al.l (|2006l ) used the non-LTE model atmo- 
sphere code CMFGEN to deduce the atmospheric parame- 
ters listed in Table [1] They estimated surface gravities at 
individual temperatures for a second order log g - T^s fit 
based on previous studies of B-type supergiants. An initial 
^ = 20 kms^^ was assumed, with values considered in the 
range of 10 - 40 kms~^ in multiples of 5 kms~^, if the 
fits of helium and silicon lines were not in agreement. Val- 
ues of the projected rotat ional velocity {v sini) were taken 
from lHowarth et al.l l|l997l) . with the exception of HD 14956 
which was not included in thei r samp le. For their temper- 
ature estimate, ICrowther et al.l (|2006l ) employed the silicon 
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Table 1 . Adopted atmospheric parameters for the sample of supergiants employed in the com parison between the model atmosphere 
calculations using C MFGEN and TLU S TY. T he spectral types and CMFGEN parameters are taken froml Crowther et alj |2006), with additional 
spectral types from iLennon et al. Note that the atmospheric parameters derived from the two codes have, in general, been 

determined using different data sets (D^). 



TLUSTY CMFGEN 



HD 


Name 


Spectral 


Teff 


log g 


5 


V sini 


Dt 


S/N 


^off 


log 9 




V sinj 


Dt 


Number 




Type 


(lO^K) 


(dex) 


(kms-i) 


(kms"-"-) 




ratio (lO^K) 


(dex) 


(kms-i) 


(kms-i) 




152235 


V900 Sco 


B0.7Ia 


22.7 


2.70 


15 


76 


1 


400 


23.0 


2.65 


10 


81 


3 


154090 


V1073 Sco 


B0.7Ia 


22.3 


2.70 


16 


80 


1 


420 


22.5 


2.65 


10 


78 


3 


148688 


V1058 Sco 


Blla 


20.6 


2.50 


16 


71 


1 


435 


22.0 


2.60 


15 


72 


3 


14956 


V475 Per 


B1.5Ia'' 


19.4 


2.40 


20 


75 


2 


370 


21.0 


2.50 


10 


80 


2 


14818 


V554 Per 


B2Ia'' 


18.8 


2.40 


18 


76 


2 


240 


18.5 


2.40 


20 


82 


2 


14143 




B2Ia'' 


19.1 


2.40 


18 


72 


2 


230 


18.0 


2.25 


20 


76 


2 


53138 


o2 CMa 


B3Ia'' 


15.2 


2.00 


18 


48 


1 


330 


15.5 


2.05 


20 


58 


4 



Dt=Data source: 1 = FEROS, MPG/ESO 2.2-m telescope; 2 = ISIS, WHT; 3 = Loral 1 K spectrograph, CTIO; 4 = 

Richardson-Brealey spectrograph, JKT. 



Table 2. Adopted atmosp heric parameters for the s ample stars (nine supergian ts and two main-sequen ce stars). Spectral types are fro m 
the Bright Star Catalogu e llHofHeit &: Warrenlll995h . with additional data from I Ferryman et al. I lll997ll °. ISchmidt fc CarruthersI lll995ll '' 
and lLennon et aU l ll992l) '^. All of the objects have been observed using FEROS on the MPG/ESO 2.2-m telescope. 



HD Name Spectral T(,ft log 5 ^ vsmi S/N 

Number Type (lO^K) (dex) (kms~i) (kms'^) ratio 



155985 




B0.5Iab'^ 


23.4 


3.00 


13 


63 


370 


108002 




Bllab" 


20.4 


2.70 


16 


65 


385 


142758 


V361 Nor 


B1.5Ia" 


19.0 


2.50 


22 


65 


415 


141318 


V360 Nor 


B2II 


20.6 


3.30 


13 


43 


475 


53138 


o2 CMa 


B3Ia'= 


15.2 


2.00 


18 


48 


330 


51309 


tCMa 


B3II 


16.6 


2.30 


19 


49 


370 


159110 




B4Ib'' 


20.1 


3.30 


10 


18 


500 


164353 


67 0ph 


B5IP 


15.8 


2.60 


23 


9 


500 


91619 


V369 Car 


B7Iae 


13.2 


1.80 


20 


29 


500 


126341 


Lup 


B2IV 


21.5 


3.65 


3 


18 


500 


79447 


t Car 


B3III 


18.0 


3.60 


9 





420 



ionization balance with lines of Sill (4128 and 4130 A), Sim 
(4552, 4568 and 4574 A) and Siiv (4089 A). The Sim and 
Siiv lines were used for BO. 7 - B2 supergiants, and Sill/Sim 
for B2.5 - B3 objects. 

The TLUSTY atmospheric parameters have been de- 
termined using an i nterre lated iterative process (see 
[Thompson et al.ll2006l . |2007| ). Effective temperatures. Toff, 
were determined using the silicon ionization balance from 
the multiplets of Sill (4128, 4130 A) and Sim (4552, 4567, 
4574 A) and the Siiv line at 41 16 A (also the 4088 A Siiv 
line for HD 126341). Similarly to lCrowther et al.1 (|2006l ). the 
Si Ill/Siiv lines were employed for stars in the spectral range 
B0.5-B1.5, whil e Siii/Sim line s were used for B2-B7 ob- 
jects. However, ICrowther et al.l (|2006l ) employed the Siiv 
line at 4088 A (which can be blended with 11 4089.29 A) 
rather than 4116 A, as used here. In the B1.5 and B2 stars, 
in addition to other silicon lines, blended or weak features of 
both Si II and Si iv were observed and the use of these lines 
is discussed below. 

Surface gravity estimates, log g, were obtained by over- 
plotting theoretical profiles on to the observed spectra in 
the region of the Balmer lines, with an associated error of 
±0.2 dex due to the uncertainty in profile fitting. Note that 



only the Balmer lines H5 and H7 have been used, as wind 
effects appe ar small for these lines but can be significant for 
Ha and H/3 (|McErlean et al.ll 19991 ). The microturbulence, 
was determined using the Si ill multiplet around ~ 4560 A, 
by plotting abundance against line strength and determining 
when the gradient is zero. 11 lines were not used as micro- 
turbulence estimates from this species yields larger values 
than t hose from Sim lines, especially in the cas e of super- 
giants (|McErlean et al.ll 19991 : [Trundle et al.ll2004h . Values of 
^ in Tables [1] and [2] have uncertainties of at least 5 km s~^. 

The effect of rotational line broadening, v smi, was cal- 
culated by overlaying instrumentally broadened theoretical 
profiles, with the appropriate atmospheric parameters, on 
to the individual lines of the Sim multiplet. Theoretical 
profiles were scaled to the same equivalent width as the 
observed features, convolved to the same spectral resolu- 
tion, and then rotationally broadened and overlain on to 
the spectra, to determine the most appropriate v sini value 
(see Hunter ct al. 2003 for further details). Other large-scale 
for ms of broadening, s uch as macroturbulence (see for e-Kam- 
ple lRyans et al.|[2002l ). were not considered due to the diffi- 
culty in distinguishing between types of broadening through 
profile fitting. Thus, the derived v sini values should be con- 
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sidered as upper limits. With the e xception of HP 53138 , 
thes e estimates agreed wit h those of iHowarth et all (|l997l ) 
and lCrowther et al.l ()2006l ) to within ± 6 kms"\ 

As stated above, in the B1.5 and B2 stars, blended or 
weak features were observed for Si ll and Si iv. In the cases 
of HD 126341, HD 14956, HD 142758 and HD 141318, the 
Sill/Silll ionization balance gave a Tea estimate approxi- 
mately 900 K larger than obtained using the Silll/Siiv lines, 
with the gravity estimate changing by ~ 0.1 dex. As all of 
the features appear reliable, intermediate values of Tcft and 
log g were considered the most appropriate for the objects 
and these are listed in Tables [T] and [2] For the other B2 
objects (HD 14818 and HD 14143), the Sin lines were used, 
as the Si iv features appeared weak and unreliable and thus 
were not employed to estimate the effective temperature. 

The atmospheric parameters of HD 152235, HD 154090, 
HD 14956, HD 14818 and HD 14143 lie at or below the lowest 
available gravity point in the metal line grid. However, using 
the method detailed above, it was possible to deduce atmo- 
spheric parameters by extrapolating the silicon abundance 
estimates, and fitting the Balmer lines using the higher res- 
olution grid. This was not believed to be a significant source 
of error, as tests performed using the available grid points in- 
dicated a difference of < 0.1 dex between silicon abundances 
determined using adjacent grid points. 

In general, good agreement is found between the at- 
mospheric parameters of the two studies, within the errors, 
despite the use of different model atmospheres. Compar- 
ing the temperature structures of CMFGEN and tlusty as a 
function of the Rosseland optical depth for HD 53138 (Fig. 
[l|, indicates that the temperature structures are in good 
agreement in the regions of formation for Felll and Sim, 
i.e. logio (Rosseland optical depth) ~ -1.8 and -1.5, respec- 
tively. However, there are differences in the temperature and 
gravity estimates for HD 148688, HD 14956 and HD 14143. 
This may be due to the use of different silicon lines, espe- 
cially as the CMFGEN and tlusty atmospheric parameters 
were deduced from the same observational data for the latter 
two objects. The v sini valu es found here tend to be lower 
than those of ICrowther et a l. (2006), with the exception of 
HD 154090. In general, the usini estimates derived here are 
thought to be more reliable due to the improved observa- 
tional data and methods of analysis. 



3.3 Iron line selection and abundances 

Felll lines were identified from synthetic spectra produced 
with the TL USTY line list, using a similar method to that 
discussed bv iGrigsbv et al.l (|l996l ). Spectra were computed 
covering the wavelength range observed with the FEROS 
data, with atmospheric parameters corresponding to a typ- 
ical B-type star of spectral type B2, e.g. HD 141318, as this 
spectral type appears to have the strongest Felll spectrum 
from our sample. (Our models indicated that Felll peaks at 
Tcft ~ 21000 K). Two sets of abundances were considered; 
one with a solar abundance for all elements, and the other 
with all elements excluded apart from iron. These were then 
overplotted on to the HD 141318 spectra and examined for 
Feiii features. This process was then repeated for the B3 
object HD 79447 to confirm the selection. This object was 
chosen despite having a weaker Fe ill spectrum than the B2 



Model At.Param 

CMFGEN 15500,2.05,20 

TLUSTY 15000,2,00,20 

--TLUSTY 15500,2.00,20 




-2 
bglO (Rosseland Optical Depth) 



Figure 1. Temperature structure plotted against Rosseland 
mean optical depth for CMFGEN and TLUSTY, with the curves rep- 
resenting the different atmospheric parameters (At.param: Tefj, 
log g and ^) for HD 53138, as calculated using the two models. 
The TLUSTY curves are for the closest available grid points to the 
derived parameters of T^g = 15200 K and g = 18 kms~^. The 
regions of formation for Fe III and Si III are indicated. 



objects, as it is a slower rotator and hence possible blended 
features could be more easily identified. 

The amount of rotational line broadening in the theo- 
retical spectra was increased to observe the effects on the 
identified lines. If two or more iron features could not be re- 
solved, these were treated as blends within tlusty. Where 
an Felll line was blended with an element of a different 
species, e.g. O ll, the iron line was included provided that 
the features could be adequately resolved from each other. 
A number of Felll lines that have been considered in pre- 
vious studies were excluded from our investigation if they 
were deemed too weak or blended with other species. In 
particular, we note the following: 

• The feature at 4395 A dete cted bvlvan HeldenI (|l972al '): 
iPintado fc AdelmanI l|l993l ) and lLiu et al.l l| 19981 ). was found 
to be blended with On at 4395.93 A. 

• A Felll line at 4168.40 A, included by 

IPintado fc AdelmanI l|l993h . is a blend of Felll (4168.45 A) 
and Sii (4168.38 A). 

• For the Fe ill lines around 4372 A it was not possible to 
resolve the individual features in the tlusty models, or in 
the slowest rotating object, HD 79447. Also, they appear to 
be blended with C ll lines. 

• There are a number of Fe ill lines in the range ~ 5880 - 
6000 A, but this spectral region contains telluric transitions, 
making identification and measurement difficult. 



As noted in Section |3T1 all of the tlusty (LTE) abun- 
dance estimates were generated using the Galactic grid and 
are given in Tables|4l[5]and[6l The CMFGEN (non-LTE) abun- 
dances (also shown in Table |4]) were obtained by calculating 
synthetic spectra for each star, using the appropriate atmo- 
spheric parameters, for iron abundances of Fe — 7.35, 7.65 
and 7.95 dex. Equivalent widths were then measured and 
compared to the observed values, and abundances obtained 
by interpolation or extrapolation. As for the tlusty models, 
this should normally not lead to significant errors. In Table 
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Table 3. Sample of the Fe in wavelengths and transitions, with atomic data taken from the Kurucz database" ( http: / /no va. astro. umd.edu| 
iHubenv 1988; Hub env fc Lanj 1995; Hubonv ct al. 1998) and Nahar & Pradhan (1996)''. Full details available on-line. 
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3 
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131 the absolute abundance estimates from the tlusty mod- 
els, for HD 152235, HD 154090, HD 14956, HD 14818 and 
HD 14143, have been extrapolated, as the atmospheric pa- 
rameters lie at or below the lowest available gravity grid 
points. However, tests using the available grid points indi- 
cate that this method yields reliable abundance estimates 
and should not lead to significant errors. 

The abundance values have associ ated errors due to 
both systematic and random errors (see lHunter et ai1l2005l 
for further details). For the tlusty comparisons, systematic 
errors arise due to the uncertainties in the atmospheric pa- 
rameters, and are estimated by varying the parameters by 
their associated errors. These error estimates have also been 
adopted for the CMFGEN comparison. Random errors are re- 
lated to the data analysis, such as the oscillator strengths, 
uncertainties in observations and line fitting errors. Thus, 
the random uncertainty is assumed to be the standard devi- 
ation of the abundance estimates from individual features, 
divided by the square root of the number of lines observed 
for that species. The total uncertainty is then taken as the 
square root of the sum of the squares of the systematic and 
random errors. 
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4 RESULTS AND DISCUSSION 

4.1 Comparison of CMFGEN and tlusty results 

As discussed in Section 13.11 there are different physical 
assumptions in the model atmosphere codes CMFGEN and 
TLUSTY. In particular, tlusty assumes plane-parallel ge- 
ometry, while CMFGEN uses spherical geometry and consid- 
ers wind properties. Additionally, non-LTE effects have been 
included in the CMFGEN Fe lll calculations. The effects of the 
different assumptions on Fe ill are considered for the sample 
of supergiants, shown in Table [T] which includes the atmo- 
spheric parameters det ermined here using tlusty and from 
ICrowther et al.l l|2006l ) . Table |4] shows the absolute abun- 
dance estimates deduced using the two model atmosphere 
codes, with the average values calculated from all of the 
observed Felll features, with the exception of those at 4053 
and 4222 A. In the CMFGEN models with higher temperature 
and lower iron abundance, these features go into emission, 
and thus no abundance estimates could be obtained for the 
BO. 7 - B1.5 objects. To maintain consistency, these lines 



Figure 2. Abundance estimates from TLUSTY plotted against 
those from CMFGEN for Felll features in all the comparison su- 
pergiants. The (■ ■ ■) lines represent the x = y (45°) relation (plus 

offsets of ± 0.2 dex), while the ( ) lines represent x = 7.5 dex 

and y = 7.5 dex (plus ± 0.2 dex offsets). 

have also been excluded in the B2 - B3 stars, but this does 
not significantly affect the average values. 

In general, the average abundance estimates are in ex- 
cellent agreement, having a typical difference of less than 
0.05 dex, with HD 152235 displaying a larger discrepancy of 
0.18 dex. Differences in the abundance estimates for the in- 
dividual lines are generally less than 0.2 dex, although there 
are some exceptions. Fig. [2] shows the abundance estimates 
from tlusty and CMFGEN for each Fe lll transition observed 
in all the supergiants. Assuming that all the stars have the 
same iron abundance, these values would be expected to be 
concentrated around a single point. To investigate possible 
sources of error two approaches have been adopted. Firstly, 
if errors due to the use of different physical assumptions pre- 
dominate, abundance estimates would be transposed either 
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horizontally or vertically. Given the simplification adopted 
in the tlusty calculations, horizontal shifts would appear 
more probable. Alternatively, if errors were present in both 
analyses (for example, due to incorrect oscillator strengths), 
points would be shifted along a 45° line. To illustrate the 
former, horizontal and vertical loci at an abundance of 7.5 
dex (typical of the Galactic value) plus offsets by ±0.2 dex 
are illustrated, while for the latter, loci at 45° plus offsets 
by ±0.2 dex are shown. 

From Fig. [2j the majority of points appear to lie along 
the 45° line. However, there are some estimates which fall 
outside of this relation. These are due to the 5833 A line (in 
HD 154090 and HD 148688), and most of the features in the 
spectral range 4419-5193 A observed in the two BO. 7 stars 
(HD 152235 and HD 154090). This suggests that the differ- 
ences in physical assumptions may become more significant 
at higher temperatures and longer wavelengths. High abun- 
dance estimates situated along the 45° relation include the 
4273 and 4296 A features, which are from the same multi- 
plet, suggesting that the oscillator strengths of this multiplet 
may be too low. Conversely, estimates due to the 4005 A line 
for almost all of the supergiants are low, strongly suggesting 
that the adopted oscillator strength of this transition is too 
high. There are also other anomalous points, due to lines at 
4164, 4166, 4352, 5193, 5306 and 5573 A. These features, 
when observed in other objects, have abundance estimates 
situated between the horizontal loci, thus making it difficult 
to confirm their behaviour. 

The Fe ionization structure included in the model at- 
mosphere codes CMFGEN and tlusty may not be correct, 
potentially leading to erroneous estimates of the relative 
strengths of Fe ll, Fe ill and Fe iv lines. Examining the Fe ion- 
ization structure in CMFGEN, in cooler stars (e.g. HD 53138) 
Fe^"*" is found to be dominant in the region of Fe ill line for- 
mation. However, in the hotter stars (e.g. HD 152235) Fe^^ 
and Fe"^"*" compete for dominance in the Felll line forma- 
tion region. Incorrect ionization structure should lead to a 
similar shift observed along the 45° line for all the Felll 
lines. However, such a trend is not observed in the sample 
of supergiants, suggesting that uncertainties in the ioniza- 
tion structure are not a significant source of error. 

Therefore, the major sources of error appear to be 
present in both analyses, as the points are distributed along 
the 45° line. The most likely causes are observational un- 
certainties, such as equivalent width measurements, due to 
blending and poor quality data (especially in the later spec- 
tral type objects), or incorrect log <;/ values. Differences be- 
tween the average abundances obtained from the two ap- 
proaches are relatively minor, as there are few points lying 
outside the 45° relation. Thus, it appears that in general 
the different physical assumptions, in particular non-LTE 
effects, are small and a LTE approximation may be a reli- 
able assumption for the Fe lll spectrum in B-type stars. 

4.2 Assessment of individual Felll lines 

Although the mean stellar abundance estimates agree, to 
within the errors, with the Galactic metallicity value, the 
reliability of the individual features varies. The majority of 
the Feiii features observed are due to single lines or blends 
of lines from the same multiplet, but some are blends with 
other multiplets (see Table |3]). The reliability of individual 



lines is discussed below in terms of equivalent width mea- 
surements and abundance estimates (see Tables [5] and |6} . 

• 3953 A: A weak feature only observed in the narrowest- 
lined stars (HD 159110, HD 126341 and HD 79447). Abun- 
dance estimates for this line are in excellent agreement with 
the average for these objects. 

• 4005, 4022, 4039 A: The 4005 A line is a relatively 
weak feature and is situated in the wings of He I 4009 A. Its 
abundance estimate is consistently lower than the average 
for the individual objects by more than 0.15 dex. However, 
the other two features (4022 and 4039 A) yield similar abun- 
dances, indicating that the values for this multiplet are con- 
sistent, although there may be problems with the oscillator 
strengths. The 4022 A feature in the supergiant HD 159110 
gives a larger abundance than the other features observed, 
by ~ 0.4 dex, suggesting it may not be reliable. 

• 4035, 4053, 4081 A: The 4035 A line observed in 
HD 79447 yields a larger abundance than the other features 
of this multiplet. This may be due to blending with Nil 
4035.10 A and On 4035.49 A, making the equivalent width 
measurement unreliable. The 4053 and 4081 A features are 
in good agreement, despite 4053 A being a blend of two 
lines. In general, their abundance estimates are lower than 
the average value obtained, and thus may not be reliable. 

• 4122 - 4166 A: The 4164.92 A feature is close to Sll 
4165.10 A. However, this does not appear to be a prob- 
lem, as the abundance estimate from the line agrees with 
that from other features in the same multiplet. The 4166 A 
line in HD 51309 provides a low abundance as it is a weak 
feature. All lines apart from 4154 A are blends but as the 
abundance values are consistent, this is not considered a 
significant source of error. 

• 4222 - 4261 A: The 4222 A feature consistently gives 
a larger abundance than the average for a given star, which 
may be due to this line being treated as a blend within 
TLUSTY and the difficulties in accurately determining the 
equivalent width, potentially leading to an overestimation 
of its value. Other features from the same multiplet are 
observed in HD 159110, HD 126341 and HD 79447. In the 
former two stars, the abundance estimates from these lines 
tend to be lower than that of 4222 A, but in HD 79447 the 
4248 A line is in agreement with 4222 A, suggesting the 
value may be reliable. The 4248 and 4261 A lines are con- 
sistent for HD 79447 but not HD 159110. However the lines 
are extremely weak in the supergiant and hence have large 
uncertainties in their equivalent width measurements. 

• ^^75 - 4310 A: These features are all components of 
the same multiplet. In general they tend to give abundance 
estimates higher than the overall average for the stars in 
which they are observed. In the two main-sequence stars, 
these lines agree well, while in the supergiants 4273 A tends 
to yield a larger value than both the average abundance 
estimate and those from the other lines of the multiplet. This 
may be due to the line being close to On 4273.10 A, and 
thus it is difficult to resolve in stars with broader features. 



Table 4. Felll equivalent width measurements (EW)* and absolute abundance estimates* for the supergiants derived using TLUSTY (T; LTE) and CMFGBN (C; non-LTE). 
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* Equivalent width measurements: 'a': accurate to better than ±10 per cent, 'b': to better than ±20 per cent, and 'c': to less than ±20 per cent. 

t Extrapolation of values due to atmospheric parameters being on/below the available TLUSTY grid point. 

* Logarithmic abundance [M/H] on the scale log[H] = 12.00 (dex). 
" Treated as blends in TLUSTY. 

t Calculated using all lines except those at 4053 and 4222 A. 

X No abundance estimate due to cmfgen model appearing to go into emission. 



Table 5. Felll absolute abundance estimates (Ab.*) and equivalent width measurements (EW)* for the supergiants. 
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7.42 


16 c 


7.38 


11 c 


7.39 


11 c 


7.46 


13 c 7.60 


5086.00" 


19 c 


7.53 


41 b 


7.58 


65 c 


7.75 


34 c 


7.65 


37 b 


7.53 


41 b 


7.61 


21 c 


7.47 


25 b 


7.68 


27 c 7.67 


5127.00" 


53 b 


7.45 


102 a 


7.46 


139 a 


7.57 


79 a 


7.50 


93 a 


7.40 


93 a 


7.47 


53 a 


7.31 


54 b 


7.35 


65 b 7.49 


5156.11 


62 b 


7.49 


111 a 


7.49 


145 a 


7.55 


74 a 


7.41 


93 a 


7.34 


90 a 


7.40 


43 a 


7.17 


46 b 


7.21 


56 b 7.36 


5193.91 






29 c 


7.22 


43 c 


7.36 


29 c 


7.42 


31 b 


7.25 


26 c 


7.21 


20 b 


7.27 


21 c 


7.35 


21 c 7.35 


5235.66 


27 c 


7.42 






18 c 


7.35 


32 c 


7.68 


30 c 


7.69 


26 c 


7.78 


18 b 


7.45 








5272.00" 










38 c 


7.55 


26 c 


7.66 


11 c 


7.21 


10 c 


7.33 


18 b 


7.43 


12 c 


7.64 




5276.00" 










36 c 


7.40 


31 b 


7.67 


19 c 


7.36 


23 c 


7.62 


21 b 


7.41 


27 c 


7.94 




5282.00" 










36 c 


7.28 


42 b 


7.64 


20 c 


7.27 


18 c 


7.37 


25 a 


7.41 


9 c 


7.26 


7 c 7.36 


5284.00" 










11 c 


7.37 


14 c 


7.60 










7 c 


7.36 








o^yy . yo 










23 c 


7.56 


19 c 


7.50 


12 c 


7.29 


12 c 


7 43 


13 c 


7 31 








5302.60 










24 c 


7.53 


22 c 


7.52 


12 c 


7.25 
















5306.00° 










15 c 


7.55 


21 c 


7.67 


9 c 


7.31 


6 c 


7.25 


11 c 


7.43 








5573.42 










32 b 


7.53 


10 c 


7.20 


8 c 


7.10 


8 c 


7.16 


8 c 


7.17 








5833.94 


48 b 


7.14 


73 a 


7.33 


134 a 


7.77 


80 a 


7.80 


66 a 


7.69 


60 a 


7.73 


46 a 


7.52 


26 b 


7.67 


24 a 7.99 


6032.00" 














12 c 


6.90 










13 c 


7.01 








6036.00" 














12 c 


7.73 










6 c 


7.48 








A\'cragc 




7,. 52 




7.51 




7.59 




7,56 




7,36 




7,48 




7,38 




7,47 


7,52 


± 








(1.11 








11.15 




(Lin 




([.17 




(I.Ki 




ll.2:i 





* Equivalent width measurements: 'a': accurate to better than ±10 per cent, 'b': to better than ±20 per cent, and 'c': to less than ±20 per cent. 

* Logarithmic abundance [M/H] on the scale log[H] = 12.00 (dex). 
" Treated as blends in tlusty. 
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• 4352 - 4431 A: None of the lines is a blend. The lines at 
4419 and 4431 A show good agreement for all stars, although 
in HD 155985 the 4431 A abundance is slightly higher, due 
to the difficulty in resolving this feature in the spectra from 
Nil 4431.81 A. The 4352 A line tends to give a lower abun- 
dance than the average, while 4382 A gives a higher value. 
However, the overall abundance value for this multiplet ap- 
pears reliable. The 4371 A line is only reliably resolved in 
the main-sequence objects and HD 159110 as it is close to a 
strong blended feature at 4372 A due to Felll and Cll. An- 
other component of this multiplet, at 4395.76 A, has been 
excluded as it is blended with On 4395.93 A (see Section 
[33J. 

• 4569 A: This feature is only observed in HD 79447 as it 
is intrinsically weak. The abundance obtained is lower than 
the average, and does not appear reliable. 

• 5063 - 5193 A: These are strong, well identified fea- 
tures. The abundance estimates from these lines tend to vary 
around the average and appear reliable. The 5086 A line is a 
blend of several transitions and gives a slightly larger abun- 
dance than the other lines (which are components of the 
same multiplet) for the supergiants, possibly due to poor 
S/N around the line. As this trend is not observed in the 
main-sequence stars, the line appears reliable. 

• 5235 - 5282, 5299, 5302 A: For the main-sequence stars 
there is good agreement between the abundance estimates, 
but the line at 5299 A yields a higher value for HD 79447, 
and 5276 A for HD 164353. The 5276 A feature is close to 
Fell 5276.00 A, which appears to affect the abundance es- 
timates obtained for HD 164353. The Fe iii fine at 5235 A 
is adjacent to Fell 5234.62 A. However, these lines are well 
resolved in all stars apart from the supergiants of spectral 
type B3. Thus for these objects (HD 51309 and HD 53138) 
the abundances may be higher due to the difficulty in mea- 
suring the individual equivalent widths. The blending of dif- 
ferent transitions in the lines does not appear to affect the 
abundance estimates obtained. 

• 5272, 5284, 5306 A: The 5272 and 5306 A features con- 
tain transitions from several multiplets, whereas 5284 A is 
a blend lines from the same multiplet. Good agreement is 
found among all the lines, although the 5284 A abundance is 
slightly lower for HD 126341 and HD 142758, while the 5272 
A value is slightly lower in HD 79447. This is probably due 
to the lines being weak in these stars, but overall these lines 
appear reliable. 

• 5575 A: This is a blend of two transitions, and is only 
observed in HD 126341. It has a small equivalent width and 
it gives a slightly lower than average abundance estimate. 

• 546O - 5573 A: The 5573 A line tends to give a lower 
than average abundance. The abundances from the other 
lines are also slightly below the average, but are only ob- 
served in the main-sequence stars. In general these are weak 
lines and thus may not be reliable. 

• 5833 A: The abundance from this feature varies around 
the average in the stars, despite being a strong, isolated, 
single transition, and this may be due to non-LTE effects. 

• 5999 A: This line is situated close to telluric features 
(see Section 1331 and is only resolved in HD 79447. However, 
this abundance estimate appears reliable as it agrees with 
the overall value for the star. 

• 6032, 6036 A: The 6032 A hue gives a consistently lower 
than average abundance, while the 6036 A abundance esti- 



Table 6. Feni equivalent width measurements (EW)* and abso- 
lute abundance estimates for the main-sequence stars. 



HD 126341 HD 79447 



Wavelength 

(A) 


EW 

(mA) 


Abundance® 


EW 
(mA) 


Abundance* 


3953.76 


4 c 


7.59 


3 c 


7.68 


4005.04 


17 b 


7.38 


14 c 


7.39 


4022.00" 


11 c 


7.39 


9 c 


7.44 


4035.43 






6 c 


7.90 


4039.16 






6 b 


7.50 


4053.00" 


12 b 


7.43 


6 b 


7.45 


4081.01 


14 b 


7.51 


6 b 


7.38 


4122.03"*' 


18 b 


7.55 


11 c 


7.73 


4122.78"'' 


15 b 




12 c 




4137.00" 


28 b 


7.61 


12 b 


7.39 


4139.00" 


29 b 


7.59 


13 b 


7.49 


4154.96 


7 c 


7.47 






4164.73"'' 


41 c 


7.70 


15 c 


7.50 


4164.92"'' 






7 c 




4166.00" 


15 a 


7.53 


8 c 


7.47 


4222.00" 


23 a 


8.14 


8 a 


7.80 


4238.62 


6 b 


7.43 


3 c 


7.56 


4248.00" 


8 c 


7.56 


7 c 


8.02 


4261.00" 


4 c 


7.57 


3 c 


7.88 


4273.00" 


8 c 


7.64 


5 b 


7.94 


4286.00" 


11 b 


7.73 


6 c 


7.85 


4296.00" 


13 a 


7.78 


7 b 


7.90 


4304.00" 


16 b 


7.87 


7 b 


7.81 


4310.00" 


17 b 


7.88 


8 b 


7.80 


4352.58 


10 b 


7.17 


10 c 


7.38 


4365.64 


5 c 


7.24 


7 c 


7.55 


4371.34 


20 c 


7.55 


12 c 


7.33 


4382.51 


11 c 


7.93 


11 b 


8.02 


4419.60 


36 a 


7.77 


33 a 


7.59 


4431.02 


24 b 


7.68 


23 a 


7.68 


4569.76 






3 c 


7.48 


5063.42 


11 b 


7.58 


10 c 


7.71 


5086.00" 


20 b 


7.70 


14 c 


7.66 


5127.00" 


52 a 


7.67 


55 b 


7.77 


5156.11 


40 a 


7.77 


32 a 


7.46 


5193.91 


19 b 


7.49 


14 c 


7.55 


5235.66 


15 b 


7.60 


5 b 


7.49 


5272 00" 


9 b 


7.52 


6 b 


7.58 


5276.00" 


18 b 


7.63 


6 b 


7.45 


5282.00" 


20 a 


7.60 


10 b 


7.59 


5284.00" 


6 c 


7.37 


5 c 


7.95 


5299.93 


13 b 


7.51 


9 b 


7.83 


5302.60 


15 b 


7.62 


6 c 


7.59 


5306.00" 


11 c 


7.58 


8 c 


7.91 


5375.00" 


4 c 


7.49 






5460.80 


9 c 


7.52 






5485.52 


8 c 


7.35 


5 c 


7.54 


5573.42 


9 b 


7.28 


2 c 


7.16 


5833.94 


36 a 


7.89 


15 a 


7.59 


5999.00" 






10 c 


7.69 


6032.00" 


12 b 


7.16 


6 c 


7.35 


6036.00" 


5 c 


7.38 


2 c 


7.71 


Average 
± 




7.57 
0.22 




7.63 
0.26 



* EW measurements: 'a': accurate to better than ±10 per cent, 
'b': to better than ±20 per cent, and 'c': to less than ±20 per 
cent. * Logarithmic abundance [M/H] on the scale log[H] = 
12.00 (dex). " Treated as blends in TLUSTY. '' Equivalent widths 
of adjacent lines were combined to produce the corresponding 
abundances shown. 
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Table 7. Summary of previously d eterm ined iron 
abundances from iReiter mann et alj lll990l ) (R90), 
N2emczura_fc_^aszvnska ^as^ |2005^ (N05), 

Gies fc Lambert ( 1992 1 (G92) and IProfEtt fc Quidevl pOOlh 
(POl). 



HD Abundance (dex) A[Fe/H] Ref. 

Number This work Previous Work 



79447 7.63 


7.62" 


+0.01 


R90 




7.66'' 


-0.03 


R90 


126341 7.57 


7.11= 


+0.46 


N05 


51309 7.36 


7.51'= 


-0.15 


G92 




7.69^* 


-0.33 


G92 




6.93*= 


+0.43 


POl 


" Using Fe III lines, 


Using Fe II lines, 


Teft = 17390 K, 


14860 K, Assuming 


; a Galactic abundance of 7.50 dex. 



mates have a larger scatter. These lines have relatively small 
equivalent widths, and do not appear reliable due to the low 
S/N of the region. 

As there are a number of Fe ill lines observed it is pos- 
sible to compare microturbulence values obtained using the 
Felll and Sim lines. In general, the values obtained are con- 
sistent, when using all of the Felll features observed. The 
multiplets of a^P^-z^P" (4352 - 4431 A) and b^D^-z'^P" 
(5063 - 5193 A) could be considered independently and are 
in good agreement with the overall result, provided well- 
observed lines are employed. 

4.3 Comparison with previous studies 

Most of our targets have been analysed previously. For 
HD 108002 and HD 159110, no previous results were found, 
while for HD 155985 and HD 142758 only project e d ro- 
tational velocities have been determined (iBalonal 1 19751 : 
iHowarth et al.lll997l '). Only three of the objects have previ- 
ous iron abundance estimates, namely HD 79447, HD 126341 
and HD 51309. These are discussed below and summarized 
in Table [71 

IReitermann et al. I (1990) analysed optical data taken 
with the CASPEC spectrograph on the 3.6-m telescope 
(ESO) for HD 79447, covering the 3900 - 4900 A wavelength 
range. They obtained LTE atmospheric parameters of Tcff 
= 17800 K and log g = 3.50 dex, and Fell and Felll iron 
abundances. Their Felll (and overall iron) abundance is in 
excellent agreement with that found here. 

In the case of HD 126341, lUE spectra spanning 1100 
- 3200 A were a nalysed with a best-fit procedure by 
iNiemczura fc Daszviis ka-Daszkicwicz (2005i). A metallicity 
was obtained, which is significantly lower than our value, 
possibly due to the use of UV spectra, as our atmospheric 
parameters agree well with theirs (Ten ~ 22700 K and log g 
= 3.74 dex). 

There ar e two previous est i mates for the iron abundance 
in HD 51309. iGies fc LambertI (|l992l ') analysed optical data 
from the McDonald 2.1-m telescope and coude spectrograph, 
identifying two Fell and two Felll lines. For the Felll lines, 
we found a similar equivalent width to their value for 5063 



A, but for the 5156 A feature their result is approximately 
50 per cent of that measured here (54 mA compared to 93 
mA). Also, they do not observe the line at 5193 A, which we 
have detected, possibly due to the use of higher S /N spectra 
in the present paper. This was tested by degrading theoret- 
ical spectra, calc ulated using the appro priate atmospheric 
parameters (from lGies fc Lambertlll992l ) and Fe abundance 
for the line (taken from this work), indicating that the fea- 
ture may have been unclear due to noise. Gies fc Lambert 
considered two different temperature estimates (17390 and 
14860 K) and associated abundances. They state that their 
data are insufficient to distinguish between these tempera- 
tures. However, both temperatur e estimates are consistent 
with that found here (16600 Kl. lProffitt fc QuiglevI (|200ll ) 
also investigated HD 51309, using lUE SWP spectra to ob- 
tain an LTE abundan ce. Their valu e is much lower than the 
abundance derived bv lGies fc Lam bert ( 1992) or that found 
here. As for HD 126341, this difference is believed to be due 
to the use of UV spectra rather than optical data, as their 
effective temperature (16570 K) and surface gravity (2.60 
dex) agree well with our estimates. 

It is worth noting that a selection of iron features have 
been consi dered in other studie s, but with no abundances de- 
termined. iLennon et all (|l993l ) observed five of our sample 
of supergiants (HD 14956, HD 14818, HD 14143, HD 53138 
and HD 164353). They note that several lines due to Felll 
are visible (over the range 3950 - 4950 A), but tha t they are 
relatively weak and hence unreliable. IVan HeldenI (|l972al lbh 
investigated HD 53138 and observed seventeen absorption 
features due to Fe ill in their optical (3500 - 8800 A) spectro- 
grams, stating that fourteen are blended with lines of other 
spectra and giving equivalent width measurements for thir- 
teen. While a number of these lines have been included here, 
several were excluded as they were considered unreliable due 
to blending with other features. For example, Felll 4041.86 
A appears to be blended with On 4041.95 A (see Section 
I3.3|l . For the lines include d here, the equivalent width mea- 
surements by Ivan HeldenI (|l972al ) tend to be larger. Those 
listed here are believed to be more reliable due to the im- 
proved quality of the spectra. 



5 CONCLUSIONS 

Due to their intrinsic weakness, Felll absorption lines have 
not been widely considered for use in chemical composition 
studies. Instead, the very rich ultraviolet spectral region has 
been favored. However, the results found in this paper sug- 
gest that the optical region can provide consistent results. 
The stars here display abundance estimates that agree with 
the Galactic metallicity, and are consistent with previous 
studies using optical spectra, where available. By contrast, 
previous determinations from ultraviolet spectra have fol- 
lowed the trends observed in other stud ies (for example 
iThompson et al.ll2007l : iDufton et al]|2007h . providing lower 
abundance estimates, in these cases by approximately 0.5 
dex. 

Although our study has concentrated on B-type stars 
found in the Milky Way, the optical Felll lines examined 
here can be applied to studies of B-type stars in other galax- 
ies, provided suitable S/N spectra are employed. For exam- 
ple, iTrundle et ai] l|2002h analysed a sample of B-type su- 
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Table 8. Recommended Felll lines for use as abundance diag- 
nostics. 



Line Spectral Type Line Spectral Type 



(A) 


Range 


(A) 


Range 


4419 


B0.5-B7 


5156 


B0.5-B7 


4431 


B1-B7 


5272 


B1.5-B5 


5063 


B1.5-B7 


5282 


B2-B7 


5086 


B0.5-B7 


5299 


B1.5-B4 


5127 


B0.5-B7 


5302 


B1.5-B4 



(B1.5) and HD 53138 (B3). However, the feature, along with 
others observed (see Tables [l] [5] and |6} may be suitable 
for diagnostic use, depending on the quality of the spectra 
used. We note that the stars included in this study have 
been selected due to being narrow lined and are typical of 
their individual luminosity classes. However, if objects with 
larger v sini values were employed, blending between Fe lines 
and other metal features may occur, thus, care should be 
taken when using objects with larger projected rotational 
velocities. 



pergiants in M 31, obtaining an Fe abundance from the Fe ill 
line at 4419 A. lRolleston et al] (|2002i ) analysed a sample of 
OB-type main-sequence stars from the LM C, finding an Fe 
abund ance for one object (PS 34-16), while iRoUeston et al.l 
l|2005 ') investigated a B-type dwarf from the SMC (AV304), 
finding agreement w ith other giant stars. More recently, 
iTrundle et al] l|2007l ) used similar methods to those detailed 
here and obtained Fe abundances for B-type stars in the 
Galaxy, LMC and SMC, using the two Fe ill hues at 4419 and 
4430 A. The results found were consistent with the present 
accepted metallicities of these systems. These studies indi- 
cate that the optical Felll lines can provide reliable abun- 
dance indicators in different galaxies. 

Our comparison of stars analysed using the model at- 
mosphere codes CMFGEN and tlusty generally shows lit- 
tle difference in the abundance estimates, indicating that 
the different physical assumptions, in particular non-LTE 
effects, are small for this species. Therefore, the results sug- 
gest that the optical Fe ill absorption line spectrum may be 
used with confidence in chemical composition studies, and 
an LTE analysis provid es reliable results. 

The atomic data of lNahar fc PradhanI (Il99( ), employed 
both in this paper and bv lCrowther et alTl|200(: ). appear to 
provide appropriate abundances, although there are some 
features, such as those at 4005 and 4273 A, whose log val- 
ues may be incorrect. Comparing the values in Table|3]shows 
that, for some features, e.g. the 4166.88, 4419.60 and 5272.90 
A lines, there are large differenc es between the atomic dat a 
from the Kurucz database and iNahar &: PradhanI (|l996f ). 
Further work is required to refine the atomic data for this 
species. 

In Table |5] we list recommended Fe ill lines which we 
believe, based on the present study, will provide reliable di- 
agnostics for the iron abundance in early B-type stars. Lines 
have been selected based on the following criteria: 

• Relatively strong, isolated feature, free from known 
blends. 

• Yields an abundance estimate within ±0.2 dex of the 
mean iron abundance for all well observed lines. 

The range of spectral types over which it is advisable 
to use the lines as an abundance diagnostic is also listed in 
TableO This spectral type range is not the same for all lines, 
generally due to the fact that some of the weaker transitions 
are only detected over a restricted span of spectral types. 
Observations for all of the lines in this Table are shown in 
Fig.[3]for HD 79447. The line at 5193 A has been included in 
the Figure, but not the Table, because its derived abundance 
is more than 0.2 dex larger than the mean value for three of 
the stars studied here, namely HD 108002 (Bl), HD 142768 
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Figure 3. Observed Fein lines in HD 79447 as listed in TablefS] including the 5193 A line. Overplotted are theoretical fits (smooth line) 
to the Fe III lines, calculated using the average abundance estimate for the star of 7.63 dex. 
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Table 3. Fein wavelengths and transitio ns, wi th atomic data taken from the Kurucz database" l |http: //nova. astro. umd.edu[ iHubenvl 
ll988l:lHubenv fc Lan z 1995; H ubenv et allll99^ and lNahar fc PradhanI lll996l ')''. 
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